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Abstract 

In this paper, the second law analysis of hydromagnetic 
conducting flow due to a stretching rotating disk with heat 
transfer is investigated using a semi-analytical/numerical 
technique termed DTM-Pade simulation. The study has 
applications in rotating magneto-hydrodynamic (MHD) 
energy generators for new space systems and also thermal 
conversion mechanisms for nuclear propulsion space 
vehicles. The momentum and energy conservation equations 
are non-dimensionalized using appropriate transformations 
leading to a set of nonlinear, coupled, ordinary differential 
equations for momentum in the radial, azimuthal and 
normal directions and a temperature distribution equation. 
Using the appropriate velocity components and temperature 
field, the entropy generation equation is obtained. The 
effects of various parameters such as magnetic interaction 
parameter, rotation parameter, Eckert number, Brinkman 
number on the entropy generation and Bejan number are 
illustrated and described. The irreversibility mechanisms of 
the entropy generation for the emerging parameters are also 
investigated and suggestions for minimizing the entropy 
generation are proposed. The DTM-Pade results are verified 
with the numerical results and found to be in excellent 
agreement. The simulations also show the feasibility of using 
magnetic rotating disk drives in novel nuclear space 
propulsion engines. 
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Introduction 

The second law of thermodynamics is employed by 
engineers to obtain the optimal design of thermal 



systems via minimizing the irreversibility and entropy 
generation, which can improve the efficiency of 
industrial systems. Entropy generation involves 
thermodynamic irreversibilities, such as characteristics 
of convective heat transfer, heat transfer across finite 
temperature gradients, viscous dissipation effects and 
magnetic field effects which arise in modern aerospace 
heat transfer processes. The fluid flow and heat 
transfer processes are intrinsically irreversible, which 
leads to an increase entropy generation and useful 
energy destruction. Entropy generation has stimulated 
significant interest in recent years in aerospace thermal 
sciences. Erbay et al. used second law analysis to 
investigate computationally the entropy generation in 
a channel with a finite volume method and SIMPLE 
algorithm. Aiboud and Saouli applied second law 
thermodynamics to analyze the entropy generation in 
magneto-viscoelastic flow over a stretching surface. 
Al-Odat et al. simulated the effect of magnetic field on 
the entropy generation due to laminar forced 
convection over a horizontal flat plate with an implicit 
finite difference technique. Aiboud and Saouli further 
studied second law thermodynamics and viscoelastic 
magnetized flow from a stretching surface with 
double-diffusive convective heat and mass transfer, 
elaborating the influence of magnetic parameter, 
Reynolds number and Prandtl number on the entropy 
generation number. Sahin described the effect of the 
variable viscosity on the entropy generation in a 
laminar fluid flow, showing that the entropy 
generation due to viscous friction became dominant 
for low heat-flux terms. Ibanez and Cuevas examined 
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entropy generation minimization of a hydromagnetic 
flow in a micro-channel. Makinde and Beg employed a 
perturbation expansion technique coupled with a 
special Hermite-Pade' approximations in the MAPLE 
program, to investigate volumetric entropy generation 
numbers, irreversibility distribution ratio and the 
Bejan number evolution in magneto-hydrodynamic 
channel flow in a plasma propulsion duct. Arikoglu et 
al. reported on the effect of slip in entropy generation 
from rotating disk in hydromagnetic flow with a 
differential transform method. Entropy generation 
analysis was applied to modeling and optimization of 
MHD induction devices by Salas et al. 

Rotating Magneto-hydrodynamic (MHD) flows 
exploit magnetic fields which can induce currents in a 
movable conductive fluid. Liquid metals, plasmas and 
electrolytes are all important examples of MHD fluids. 
MHD flows also utilize a Lorentzian drag force which 
can be used to regulate a variety of flow regimes. 
Many applications exist for rotating hydromagnetic 
flows including propulsion systems, rotating MHD 
energy generators, smart spacecraft landing gear 
systems, hydrogen production with solar MHD plants, 
plasma fusion technology, nuclear thermal control 
systems, magneto-hydrodynamic chemical reactor 
processing and biomagnetic reactors. Many further 
applications of rotating hydromagnetic flows 
including swirling disk flows, flows from revolving 
cones in liquid metal stirring etc, are documented in 
the recent monograph by Beg et al. 

In order to describe many physical systems with a 
mathematical model, nonlinear equations are used. 
Merical methods are frequently deployed to solve 
these complex systems of differential equations, which 
may be multi-degree and often strongly coupled. 
Although many powerful numerical techniques exist 
including finite element methods, network simulation, 
finite difference methods and CFD codes such as 
FLUENT, computational expense is always a problem. 
Therefore semi-numerical/analytical methods, such as 
the homotopy perturbation method (HPM), 
differential transform method (DTM) and homotopy 
analysis method (HAM) have found increasing 
popularity amongst aerospace fluid dynamics 
researchers. These techniques can successfully tackle 
nonlinear differential equations for a variety of 
boundary conditions, and therefore they hold 
significant promise in nuclear engineering sciences. 

In this article, the second law analysis is employed to 
study the MHD fluid flow with heat transfer due to a 



stretching rotating disk. This problem was firstly 
studied by Turkyilmazoglu numerically, without 
considering entropy generation analysis. In this paper, 
these nonlinear swirling flow and heat conservation 
equations are solved via DTM-Pade simulation. As the 
second law analysis is more reliable than the first law 
analysis for many aerospace and nuclear engineering 
systems, entropy generation analysis is also included. 
In the current study, the effect of various 
thermophysical parameters, such as magnetic 
interaction parameter, rotation parameter, Eckert 
number, Prandtl number, Brinkman number and 
Reynolds number on velocity and temperature fields 
and also on the entropy generation and Bejan number 
are investigated in detail. The study has important 
applications in thermal optimization of novel MHD 
coupled nuclear space propulsion systems using 
rotating disks and hydromagnetics. 

Mathematical Transport Model 

The 3-dimensional steady magneto-hydrodynamic 
laminar boundary-layer flow of a Newtonian, 
electrically-conducting, viscous fluid from a rotating 
disk in the presence of an externally applied axially- 
directed uniform magnetic field is considered. This 
regime is often also termed hydromagnetic Von 
Karman swirling flow. The governing equations; the 
continuity, Navier-Stokes (momentum conservation), 
generalized Ohm's law and energy conservation 
equations can be presented, in the presence of Ohmic 
(Joule) heating and viscous dissipation, in vectorial 
form, as follows: 

V.m = (1) 
p(u-V)u = -VP + juV 2 u+J xB (2) 
J=a[E + uxB] (3) 
pc p (u-V)T = kV 2 T + >u<S>+J 2 /a (4) 

where all parameters are defined in the nomenclature. 
The disc radius is much larger than the boundary layer 
thickness, so that edge effects can be neglected. The 
applied magnetic field is steady and sufficiently weak 
to neglect the induced magnetic field i.e. a low 
magnetic Reynolds number is assumed. This is a 
reasonable assumption for the flow of certain working 
fluids in novel nuclear-MHD propulsion systems and 
also in nuclear engineering liquid metals, e.g., liquid 
sodium. The physical regime is depicted in Fig. 1 with 
reference to non-rotating cylindrical polar coordinates 
where all parameters are defined in the nomenclature. 
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The disc radius is much larger than the boundary layer 
thickness, so that edge effects can be neglected. The 
applied magnetic field is steady and sufficiently weak 
to neglect the induced magnetic field i.e. a low 
magnetic Reynolds number is assumed. This is a 
reasonable assumption for the flow of certain working 
fluids in novel nuclear-MHD propulsion systems and 
also in nuclear engineering liquid metals, e.g., liquid 
sodium. The physical regime is depicted in Fig. 1 with 
reference to non-rotating cylindrical polar coordinates 
[r,e,z). 
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FIG. 1 3-D SWIRLING VON KARMAN MAGNETO- 
HYDRODYNAMIC FLOW REGIME AND COORDINATE SYSTEM 

The disk rotates with a constant angular velocity (q) 
and also stretches with a constant rate (s) in the radial 
direction. The external uniform magnetic field B 
applies transverse to the disk plane i.e. in the z- 
direction and possesses a constant magnetic flux 
density Bo. Furthermore all thermophysical and fluid 
properties are assumed to be constant and radial and 
tangential electrical currents are neglected. The disk 
itself is electrically non-conducting, so that current 
density vanishes both at the disk surface and in the 
fluid regime. The penultimate and final terms on the 
right hand side of eqn. (9) designate the viscous 
heating and Joule heating contributions, respectively, 
and have been shown to be significant in propulsion 
applications. Following Turkyilmazoglu, the non- 
dimensional form for the mean flow velocities and 
temperature distribution are provided by Von 
Karman's classical transformations. With the aid of a 
dimensionless normal distance from the wall, 
rj = zR, and an appropriate Reynolds number 

incorporating the disk stretching rate, R = (s/v)^ 2 we 
have: 

{u^wJ^rsF^rsGi^H^),^ +{T w -T^)! (5) 

where the dimensionless functions F , G , H , and 
satisfy the following ordinary differential equations 



defining the radial, azimuthal, axial momentum 
conservation and energy conservation in the regime: 



H' + IF = 



F" -F 2 +G 2 -HF' -MF 



G"-2FG-HG'-MG = 

0" - Pr H0' + Ec[f' 2 +G' 2 )+ EcM (f 2 + G 2 ) = 



(6) 
(7) 
(8) 
(9) 



The boundary conditions at the stretching disk surface 
and in the free stream, are prescribed respectively, as: 



H = F-l = G-o) = 0-1 = 0, 
F = G = = as rj — > oo, 



at 



7 = 0, 



(10) 



Here co = Cl/s denotes the "rotation strength 
parameter" expressing the ratio of disk swirl to disk 
radial stretch. The case of co = corresponds to pure 
stretching of the disk without rotation. The nonlinear, 
coupled ODEs defined in eqns. (6)-(9) under 
boundary conditions (10) constitute a robust two-point 
boundary value problem, which may be readily solved 
with a variety of numerical or semi-numerical 
techniques. In the present study we implement a 
combination of the differential transform method 
(DTM) and Pade approximants, namely DTM-Pade 
simulation, which is described in due course. 

Entropy Generation Analysis 

The volumetric rate of local entropy generation, for the 
present hydromagnetic swirling heat transfer problem, 
in the presence of the axial symmetry, can be 
presented as follows: 



r> m 

gen 



where: 



[vr] 2 = 



-[vr] 2 +^-<s +— [(j-qv)-(e+Vxb)] (li) 
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J = <j[E + u x B] 



(13) 



(14) 



It is assumed that electric force per unit charge is 
negligible compared with V*B in Eqs. (11) and (14) 
and furthermore that electric current is much greater 
than QV. With these assumptions, Eq. (11) can be 



31 



www.fae-journal.org 



Frontiers in Aerospace Engineering (FAE) Volume 2 Issue \, February 2013 



shown to take the form: 
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The right hand side of the above entropy generation 
equation consists of three parts. The first term is the 
local entropy generation due to heat transfer 
irreversibility. The second larger part refers to the 
fluid friction irreversibility. The final part denotes 
magnetic field effects. The dimensionless form of the 
entropy generation rate, termed entropy generation 
number, defines the ratio between the actual entropy 
generation rate 5™ and the characteristic entropy 

generation rate Sq. By using the Von Karman 
transformation parameters, as defined in Eq. (5), the 
entropy generation number (n g ) becomes, for the 
present problem: 

2 \Br _f v 2Ec „,f v 



N c =a9'(rjf + ^F(rj) 



H [t]) 



Re Re " (16) 

+ Br\F'(rif + G'(t]f)+ BrM\F(rff + G(t]f) 

Inspection of Eq. (16), indicates that with increasing 
value of a, the effect of heat transfer irreversibility 
increases. As the Brinkman number and Eckert 
number increase or the Reynolds number decreases, 
the entropy generation due to fluid friction 
irreversibility increases. Irreversibility mechanism 
domination is an important aspect of entropy 
generation analysis since the entropy generation 
number does not provide any information regarding 
this. The Bejan number and the irreversibility 
distribution ratio are introduced to overcome this 
shortcoming. The Bejan number (Be) embodies the 
ratio of entropy generation which is caused through 
heat transfer to the total entropy generation. The 
irreversibility distribution ratio (cp) expresses the ratio 
between entropy generation due to fluid friction and 
joule dissipation in heat transfer. These two 
parameters, for the present swirling 
magnetohydrodynamic thermal problem can be 
expressed as: 



a ff(„f + ^F( n f + ^fw'M 2 + Sr[ Fin? + G'(,PJ + Mf f pfgf + of,) 2 



(17) 



(18) 



As illustrated by the definitions of the above 
dimensionless parameters, the behavior of the Bejan 
number and irreversibility distribution ratio are 
almost the same. Thus we only examine the Bejan 
number effects in this article. The Bejan number falls 
in the range, < Be < 1 . When Be = , fluid friction 
irreversibility dominates. For Be = 0.5 , heat transfer 
and fluid friction irreversibility effects are the same. 
For Be = 1 , the irreversibility mechanism is dominated 
by heat transfer effects. 

DTM- Pade Simulation 

In the current article, DTM is employed in conjunction 
with Pade approximants to solve the governing 
equations (6) to (9), under boundary conditions (10). 
Zhou was the first one who developed DTM for 
electrical circuit modeling. The principal attraction of 
this technique is that it can be applied directly to 
nonlinear differential equations without requiring 
discretization and linearization, a strong advantage 
over purely numerical methods such as finite elements 
and finite differences. DTM also does not require any 
perturbation parameters and this offers a great 
advantage over asymptotic and perturbation 
expansion techniques. This method has been 
successfully implemented in a diverse array of 
including nanofluid, hypersonic, swirl flow and 
geothermic flow. In order to solve highly nonlinear 
differential equation systems, arising frequently in 
nuclear engineering fluid dynamics, it is pertinent to 
combine DTM with Pade approximations and use 
DTM-Pade simulation. This method effectively 
increases the convergence of DTM. Further details are 
given in ref. Numerous mathematical techniques exist 
to increase the convergence radius of a given series. 
The Pade approximant is a rational fraction. Details of 




FIG. 2 RADIAL VELOCITY COMPONENT, F(rj) OBTAINED BY 
DTM FOR DIFFERENT N AND ORDERS OF DTM-PADE IN 
COMPARISON WITH REF. [26] FOR m = \,M =2. 
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such approximants are provided in Baker and Graves- 
Morris. The approximants are generated numerically 
using MATHEMATICA software. Stability, consistency 
and convergence are successfully achieved by DTM- 
Pade simulation. DTM alone diverges (figure 2). 

Results and Discussion 

To verify the validity of the present computations, 
results have been benchmarked with the non-entropy 
study of Turkyilmazoglu. Excellent correlation 
between the semi-numerical/analytical results 
obtained by DTM-Pade (Pade approximant of order) 
and the Chebsyhev spectral collocation numerical 
computations of Turkyilmazoglu is achieved, as 
observed in Fig 3. Confidence in the present DTM- 
Pade solutions is therefore high. We note that the case 
of equal rotational strength and disk stretching 
corresponds to co = 1. The cases of a> = i.e. pure 
stretching of the disk without rotation, and <»— »<x> i.e. 
pure rotation of the disk without stretching are 
extreme scenarios and not considered in the present 
study. 




FIG. 3 RADIAL, AZIMUTHAL AND AXIAL VELOCITY 
FUNCTIONS AND TEMPERATURE FUNCTION OBTAINED BY 
DTM-PADE (PADE APPROXIMANTS [10, 10]) IN COMPARISON 
WITH REF. [26] WITH Pr = Ec = 1, M = 2 AND co = 3 

Figures 4-7 present the velocity contours in all 
directions. As the radial coordinate increases, the 
primitive radial (u) and azimuthal (v) velocity 
components clearly increase. These velocity 
components are maximized near the surface of the 
disk i.e. at low values of axial coordinate (bottom right 
hand corner of both Figs. 4 and 5). The velocity vectors 
are shown in Fig. 7, in order to have a better grasp of 
the fluid flow. In both cases although co has been 
prescribed as 1, the disk stretch rate is very high at 10. 
Conversely in Fig. 6 we observe that axial velocity 
component is maximized near the disk surface for all 
values of radial coordinate (lower red band). The axial 
velocity clearly decays as we depart from the surface 
and is minimized at greater distances from the disk 
(z ~ 0.4) . Figure 7 distinctly shows the convergence of 



velocity vectors (u, w) towards the bottom right hand 
corner of the plot i.e. fluid is clearly drawn in a fan like 
mechanism outwards along the radial coordinate and 
in the negative axial direction. This pattern typifies 
Von Karman swirling flow. 




FIG. 4 RADIAL VELOCITY CONTOURS WITH M = m = 1, 

Pr = Ec = 2, v = 9 x 10~ 6 AND s = 10. 




FIG. 5 AZIMUTHAL VELOCITY CONTOURS WITH M = m = 1, 
Pr = Ec = 2,v = 9xl0~ 6 AND S = 10. 




FIG. 5 AXIAL VELOCITY CONTOURS WITH M = a = 1, 
Pr = Ec = 2, v = 9 x 10~ 6 AND * = 10. 




FIG. 7 VECTOR VARIABLES OF u AND w WITH M = m = 1, 
Pr = £c = 2, v = 9xl(T 6 AND * = 10. 

Figures 8-13 illustrate the influence of the emerging 
thermophysical parameters on entropy generation 
number, and Bejan number distributions, with radial 
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coordinate. In each of these graphs, dimensionless 
temperature difference is set to unity. 

The entropy generation number and Bejan number 
variation with changing magnetic interaction 
parameter are presented in Figs. 8 and 9. As the 
magnetic interaction parameter increases, the amount 
of entropy generation increases considerably (it is 
minimized in the absence of a magnetic field, M =0 ). 
This implies that in order to control the entropy which 
is generated from swirling flow on the disk surface, 
the value of magnetic interaction parameter should be 
reduced, an issue of interest in nuclear-MHD rotating 
disk propulsion. The fluid friction irreversibility is 
dominated mechanism of irreversibility, and is distant 
from of the disk surface for small values of magnetic 
interaction parameter. As the magnetic interaction 
parameter is elevated, the heat transfer mechanism 
dominates irreversibility. It should also be mentioned 
that a minimum point arises for Bejan number near the 
disk surface for any value of magnetic interaction 
parameter. At large values of M the Bejan number 
escalates and after ascending very sharply from the 
disk surface converges asymptotically to the 
maximum allowable value of unity, a trend sustained 
into the free stream. In essence when entropy 
generation is maximized, the Bejan number is 
minimized and vice versa in the swirling flow regime. 




FIG. 8 THE VARIATION OF ENTROPY GENERATION NUMBER 
N c WITH MAGNETIC INTERACTION PARAMETER FOR 
Pr = Ec = 2, Br = 4, Re = 8 AND co = 1. 




Fig. 9 THE VARIATION OF BEJAN NUMBER, Be WITH 
MAGNETIC INTERACTION PARAMETER FOR 
Pr = Ec = 2, Br = 4, Re = 8 AND to = 1. 

Figures 10, 11 illustrate the effect of rotation strength 



parameter on the entropy generation and Bejan 
number distributions with radial coordinate. An 
increase in co generates a similar but more 
pronounced response to that of increasing magnetic 
interaction parameter i.e. it causes a considerable 
increase in N G near the disk surface, leading to a 
maximum entropy generation value always at the disk 
surface itself. Further from the disk N c profiles decay 
sharply to vanish in the free stream. The Bejan number 
is generally minimized closer to the disk and 
maximized further from the disk, in direct contrast to 
the entropy generation number. In all cases, the lowest 
value of Bejan number arises in the same vicinity as 
the maximum value of entropy generation number i.e. 
near the disk surface. For larger values of rotation 
strength parameter Bejan number achieves the 
asymptotic profile sooner. In the asymptotic state 
where the profiles all converge on unity in the free 
stream, the heat transfer irreversibility is dominant. 




FIG. 10 THE VARIATION OF ENTROPY GENERATION NUMBER, 
N c WITH ROTATION STRENGTH PARAMETER FOR 

Pr = Ec = 1, M = 2, Br = 3 AND Re = 6. 



<3S 




FIG. 11 THE VARIATION OF BEJAN NUMBER, Be WITH 
ROTATION STRENGTH PARAMETER FOR 

Pr = Ec = 1, M = 2, Br = 3 AND Re = 6. 

Figures 12 and 13 depict the response of the entropy 
numbers to Eckert number. Increasing Eckert number 
significantly increases the entropy generation number, 
N c increases, as the Eckert number increases. It 
should also be noted that at a certain distance from the 
disk surface (77 = 0.45), the amount of generated 
entropy converges on the same value for all Eckert 
numbers. The influence of Eckert number on the Bejan 
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number (Be) profiles is similar to those induced by the 
rotational strength parameter. Brinkman number, 
Br = /j s 2 r 2 /kAT represents the ratio of direct heat 
conduction from the disk surface to the viscous heat 
generated by shear in the boundary layer. It is distinct 
from Eckert number which symbolizes the kinetic 
energy of flow to the boundary layer enthalpy 
difference. It is in fact more appropriate for industrial 
flows, whereas Eckert number often arises in high- 
speed aerodynamic flows. Both parameters are 
quantifications of viscous dissipation effects. Bejan 
number decays from a high value at the disk surface to 
a minimum near the disk surface and then ascends 
rapidly to converge on unity in the free stream. A rise 
in Ec clearly enhances Bejan number. In all cases we 
have considered positive Ec values i.e. heat is 
transferred from the disk surface to the fluid by 
convection currents. Further from the disk surface 
heat transfer irreversibility dominates. For large 
values of Prandtl number, the fluid friction 
irreversibility dominates compared with heat transfer. 
Bejan number is clearly a maximum when entropy 
generation number is a minimum. 




n 



FIG. 12 THE VARIATION OF ENTROPY GENERATION NUMBER, 
N c WITH ECKERT (VISCOUS HEATING) NUMBER FOR 
Pr = 1, (o = M = 2, Br = 4 AND Re = 6. 




FIG. 13 THE VARIATION OF BEJAN NUMBER, Be WITH 
ECKERT (VISCOUS HEATING) NUMBER FOR 
Pr = 1, (o = M = 2, Br = 4 AND Re = 6. 

We further note that the correct trends for entropy 
generation number and Bejan number have been 
achieved by DTM-Pade simulation and have been 
concurred with the observations documented in other 



studies, for example Aiboud and Saouli and Arikoglu 
et al., the latter considering Von Karman swirling 
hydromagnetic flow with slip effects. As the Brinkman 
number increases, the effect of fluid friction 
irreversibility is amplified on the disk surface. Further 
from the disk surface, the effect of heat transfer 
irreversibility becomes larger as the Brinkman number 
decreases. For high Br , the effects of heat transfer and 
fluid friction irreversibilities are almost equal. 
Effectively, the present semi-numerical computations 
testify that the contribution of Brinkman number, as 
originally proposed and elucidated by Bejan and 
further implemented in hydromagnetic heat transfer 
by Makinde and Beg, is therefore of significance in 
second law thermodynamic analysis, and plays an 
important potential role in nuclear MHD spacecraft 
engine thermodynamic optimization. 

Conclusions 

In the current investigation, second law 
thermodynamic analysis has been utilized to obtain 
the entropy generation equations for swirling Von 
Karman hydromagnetic flow with heat transfer from a 
stretching rotating disk. A semi-numerical technique 
amalgamating the differential transform method 
(DTM) and Pade approximants i.e. DTM-Pade 
simulation, has also been applied to determine 
solutions for the transformed momentum and energy 
conservation equations describing the swirling fluid 
dynamics, under appropriate boundary conditions. 
The nonlinear coupled multi-degree boundary value 
problem has been very efficiently solved with 
benchmarking to numerical shooting quadrature, 
demonstrating excellent correlation between both 
methods. The velocity contours in all directions and 
velocity vectors have been visualized to illustrate the 
actual dynamics of the flow in primitive coordinates. 
Furthermore the influence of Brinkman number on 
entropy generation number and Bejan number, have 
also been presented and analyzed. The powerful 
contributions of disk swirl and disk stretch 
(investigated via the rotational strength parameter), 
magnetic field retardation and viscous heating on the 
flow variables have clearly been identified in the 
computations. The fundamental objective of second 
law thermodynamics analysis has been also achieved, 
namely, minimization of entropy in the swirling disk 
flow regime. The present computations have provided 
some further insight into the thermodynamics and 
fluid mechanics of proposed rotating disk MHD 
systems coupled with nuclear space propulsion 
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engines, although these require further evaluation, in 
particular for transient and magnetic induction effects, 
aspects which are under consideration by the authors. 
DTM-Pade simulation holds excellent potential in 
analyzing such aerospace problems. 
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